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ABSTRACT 

We have proposed that the first phase of steUar evohition in the history of the Universe may be 
Dark Stars (DS), powered by dark matter heating rather than by nuclear fusion, and in this paper 
we examine the history of these DS. The power source is annihilation of Weakly Interacting Massive 
Particles (WIMPs) which are their own antiparticles. These WIMPs are the best motivated dark 
matter (DM) candidates and may be discovered by ongoing direct or indirect detection searches (e.g. 
FERMI/GLAST) or at the Large Hadron Colhder at CERN. A new steUar phase results, powered by 
DM annihilation as long as there is DM fuel, from millions to billions of years. We build up the dark 
stars from the time DM heating becomes the dominant power source, accreting more and more matter 
onto them. We have included many new effects in the current study, including a variety of particle 
masses and accretion rates, nuclear burning, feedback mechanisms, and possible repopulation of DM 
density due to capture. Remarkably, we find that in all these cases, we obtain the same result: the first 
stars are very large, 500-1000 times as massive as the Sun; as well as puffy (radii 1-10 A.U.), bright 
(10^ — IO^Lq), and cool (Tgurf <10,000 K) during the accretion. These results differ markedly from 
the standard picture in the absence of DM heating, in which the maximum mass is about 140 Mq 
(McKee & Tan 2008) and the temperatures are much hotter {Tgurf >50,000K). Hence DS should be 
observationally distinct from standard Pop III stars. In addition, DS avoid the (unobserved) element 
enrichment produced by the standard first stars. Once the dark matter fuel is exhausted, the DS 
becomes a heavy main sequence star; these stars eventually collapse to form massive black holes that 
may provide seeds for the supermassive black holes observed at early times as well as explanations for 
recent ARCADE data (Seiffert et al. 2009) and for intermediate mass black holes. 
Subject headings: Dark Matter, Star Formation, Accretion 



1. INTRODUCTION 

Spolyar et al. (2008; hereafter Paper I) proposed a new 
phase of stellar evolution: the first stars to form in the 
Universe may be Dark Stars, powered by dark matter 
heating rather than by nuclear fusion. Here dark mat- 
ter (DM), while constituting a negligible fraction of the 
star's mass, provides the energy source that powers the 
star. The first stars in the Universe mark the end of the 
cosmic dark ages, provide the enriched gas required for 
later stellar generations, contribute to reionization, and 
may be precursors to black holes that coalesce and power 
bright early quasars. One of the outstanding problems in 
astrophysics is to investigate the mass and properties of 
these first stars. Depending on their initial masses, the 
stars lead very different lives, produce very different ele- 
mental enrichment, and end up as very different objects. 
In this paper, we find the mass, luminosity, temperature, 
and radius of dark stars as a function of time: our re- 
sults differ in important ways from the standard picture 
of first stars without DM heating. 

The DM particles considered here are Weakly Inter- 
acting Massive Particles (WIMPS), such as the lightest 
supersymmetric particles or Kaluza-Klein dark matter 
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(for reviews, see Jungman et al. 1996; Bertone et al. 
2005; Hooper & Profumo 2007). The search for these 
particles is one of the major motivations for the Large 
Hadron Collider at CERN. The particles in consideration 
are their own antiparticles; thus, they annihilate among 
themselves in the early universe and naturally provide 
the correct relic density today to explain the dark mat- 
ter of the universe. The first stars are particularly good 
sites for DM annihilation because they form in a very 
high density DM environment: they form at early times 
at high redshifts (density scales as (1 -I- z)^) and in the 
high density centers of DM haloes. The first stars form 
at redshifts z ~ 10 - 50 in DM halos of IO^Mq (for re- 
views see e.g. Ripamonti & Abel 2005; Barkana & Loeb 
2001; Yoshida et al. 2003; Bromm & Larson 2004; see 
also Yoshida et al. 2006) One star is thought to form 
inside one such DM halo. These halos consist of 85% 
DM and 15% baryons in the form of metal-free gas made 
of H and He. Theoretical calculations indicate that the 
baryonic matter cools and collapses via H2 cooling (Pee- 
bles & Dicke 1968; Matsuda et al. 1971; HoUenbach & 
McKee 1979) into a single small protostar (Omukai & 
Nishi 1998) at the center of the halo. 

As our canonical values, we use the standard {av) = 
3 X 10~^^cm'^/s for the annihilation cross section and 
= 100 GeV for the WIMP particle mass. In this 
paper we examine a variety of possible WIMP masses 
and cross sections, and find that our main result holds 
independent of these properties. We consider — 1 
GeV, 10 GeV, 100 GeV, 1 TeV, and 10 TeV. This variety 
of WIMP masses may be traded for a variety of values 
of the annihilation cross section, since DM heating scales 
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as {av) /m^. 

Paper I found that DM annihilation provides a power- 
ful heat source in the first stars, a source so intense that 
its heating overwhelms all cooling mechanisms; subse- 
quent work has found that the heating can dominate over 
fusion as well once it becomes important at later stages 
(see below). Paper I outlined the three key ingredients 
for Dark Stars: 1) high dark matter densities, 2) the an- 
nihilation products get trapped inside the star, and 3) 
DM heating wins over other cooling or heating mecha- 
nisms. These same ingredients are required throughout 
the evolution of the dark stars, whether during the pro- 
tostellar phase or during the main sequence phase. 

WIMP annihilation produces energy at a rate per unit 
volume 

Qdm = {<Jv)p\/m^, (1) 

where is the mass density of the WIMPs. Once the 
gas density of the collapsing protostar exceeds a critical 
value, most of the annihilation energy is trapped in the 
star. For a 100 GeV particle, when the hydrogen density 
reaches ^ 10^^ cm~^, typically 1/3 of the energy is lost to 
neutrinos that escape the star, while the other 2/3 of the 
energy is trapped inside the star. Hence the luminosity 
from the DM heating is 

Ldm ^^j QoMdV (2) 

where dV is the volume element. From this point for- 
wards, the DM heating beats any cooling mechanisms. 
This point is the beginning of the life of the dark star, a 
DM powered star which lasts until the DM fuel runs out. 
The DM constitutes a negligible fraction of the mass of 
the dark star (less than one percent of the total stellar 
mass) and yet powers the star. This " power of darkness" 
is due to the 67% efficiency of converting the DM mass 
into an energy source for the star, in contrast to the ^1% 
efficiency of converting nuclear mass into energy in the 
case of fusion. 

In this paper we ascertain the properties of the dark 
star as it grows from an initially small mass at its incep- 
tion to its final fate as a very large star. At the beginning, 
when the DM heating first dominates inside the star, the 
dark star mass for our canonical case of 100 GeV mass 
particles is 0.6Mq with a radius of 17 AU (Abel et al. 
2002; Gao et al. 2007). Then further matter can rain 
down upon the star, causing the dark star to grow. By 
the time the star reaches about 3 Mq , with a radius of 1- 
10 AU, it can be approximated as an object in thermal 
and hydrostatic equilibrium, with dark matter heating 
balancing the radiated luminosity. There is ~ lOOOM© 
of unstable gas (^ Jeans mass) which can in principle ac- 
crete onto the star, unless there is some feedback to stop 
it. We build up the dark star from a few solar masses by 
accretion, in steps of one solar mass. At each stage we 
find the stellar structure as a polytrope in hydrostatic 
and thermal equilibrium. We continue this process of 
building up the DS until the DM runs out. 

As seen in equation (1), a key element in the anni- 
hilation is knowledge of the DM density inside the DS. 
The initial DM density profile in the W^Mq halo is not 
enough to drive DM heating. We consider three mecha- 
nisms for enhancing the density. First, as originally pro- 
posed in Paper 1 and confirmed in Freese et al (2008b), 



adiabatic contraction (AC) drives up the density. As the 
gas collapses into the star, DM is gravitationally pulled 
along with it. Second, as the DS accretes more baryonic 
matter, again DM falls in along with the gas. Eventually 
the supply of DM due to this continued AC runs out, and 
the star may heat up until fusion sets in. However, there 
is the possibility of repopulating the DM in the DS due to 
a third source of DM, capture from the ambient medium 
(locco 2008; Freese et al. 2008c). Capture only becomes 
important once the DS is already large (hundreds of solar 
masses), and only with the additional particle physics in- 
gredient of a significant WIMP / nucleon elastic scattering 
cross section at or near the current experimental bounds. 
Capture can be important as long as the DS resides in 
a sufficiently high density environment. In this case the 
DS lives until it leaves the womb at the untouched cen- 
ter of a lO'^Afo DM halo. Unfortimately for the DS, 
such a comfortable home is not likely to persist for more 
than 10-100 Myr (if the simulations are right), though in 
the extreme case one may hope that they survive to the 
present day. In this paper we consider all three mecha- 
nisms for DM enhancement, including DM capture. We 
remark at the outset that AC is absolutely essential to 
the existence of the DS. None of these three mechanisms 
would operate were it not for some AC driving up the 
density in the vicinity of the DS. 

Previously (Freese et al. 2008a), we considered one 
specific case, that of 100 GeV particle mass and a con- 
stant accretion rate of 2 x 10~"^M©/yr with no capture, 
assuming that the DS stellar structure can be described 
by a series of n = 3/2 polytropes. We found that it 
takes ^ 0.5 million years to build up the star to its fi- 
nal mass, which we found to be very large, ^ SOOMq. 
While the accretion was in process, the dark stars were 
quite bright and cool, with luminosities L ~ IO^Lq and 
surface temperatures Tgurf < 10,000/r. These results 
are in contrast to the standard properties of the Popula- 
tion III stars (the first stars). Standard Pop III stars are 
predicted to form by accretion onto a much smaller pro- 
tostar (1O~^M0; Omukai & Nishi 1998), to reach much 
hotter surface temperatures Tg^^f »50,000 K during ac- 
cretion, and at the end of accretion to be far less massive 
(^ 100 — 200Mq) due to a variety of feedback mech- 
anisms that stop accretion of the gas (McKee and Tan 
2008). These differences in properties between dark stars 
and standard Pop HI stars lead to a universe that can 
look quite different, as discussed below. 

In this paper we present the results of a complete study 
of building up the dark star mass and finding the stel- 
lar structure at each step in mass accretion. We include 
many new effects in this paper compared to Freese et al. 
(2008a). (1) We generalize the particle physics parame- 
ters by considering a variety of particle masses or equiv- 
alently annihilation cross sections. (2) We study two 
possible accretion rates: a variable accretion rate from 
Tan & McKee (2004) and an alternate variable accretion 
rate found by O'Shea & Norman (2007). (3) Whereas 
previously we considered only n = 3/2 polytropes, which 
should be used only when convection dominates, we now 
allow variable polytropic index to include the transition 
to radiative transport (n = 3). (4) We include gravity 
as a heat source. In the later stages of accretion, as the 
DM begins to run out and the star begins to contract, 
the gravitational potential energy is converted to an im- 
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portant power source. (5) We include nuclear burning, 
which becomes important once the DM starts to run out. 
(6) We include feedback mechanisms which can prevent 
further accretion. Once the stellar surface becomes hot 
enough, again when the DM is running out, the radia- 
tion can prevent accretion. (7) We include the effects 
of repopulating the DM inside the star via capture from 
the ambient medium. The result of this paper is that, 
for all the variety of parameter choices we consider, and 
with all the physical effects included, with and without 
capture, our basic result is the same. The final stellar 
mass is driven to be very high and, while DM reigns, the 
star remains bright but cool. 

Other work on dark-matter contraction, capture, and 
annihilation in the first stars (Freese et al. 2008c; Taoso 
et al. 2008; Yoon et al. 2008; locco et al. 2008; Ri- 
pamonti et al. 2009) has concentrated on the case of 
the pre-main-sequence, main-sequence, and post-main- 
sequence evolution of stars at fixed mass, in contrast to 
the present work, which allows for accretion from the 
dark-matter halo. 

We also cite previous work on DM annihilation in to- 
day's stars (less powerful than in the first stars): Krauss 
et al. (1985); Bouquet & Salati (1989); Salati & Silk 
(1989); Moskalenko & Wai (2007); Scott et al. (2007); 
Bertone & Fairbairn (2007); Scott et al. (2009). 

2. EQUILIBRIUM STRUCTURE 

2.1. Basic Equations 

At a given stage of accretion, the dark star has a pre- 
scribed mass. We make the assumption that it can be 
described as a polytrope in hydrostatic and thermal equi- 
librium. The hydrostatic equilibrium is defined by 

dP GMr dMr , 2 . ^ .o^ 

and by the polytropic assumption 

P = Kp^+^/'' (4) 

where P is the pressure, p is the density, Mr is the 
mass enclosed within radius r, and the constant K is 
determined once the total mass and radius are specified 
(Chandrasekhar 1939). Pre-main-sequence stellar mod- 
els are adequately described by polytropes in the range 
n = 1.5 (fully convective) to n = 3 (fully radiative). 
Given P and p at a point, the temperature T is defined 
by the equation of state of a mixture of gas and radiation 

P(.) = ^5^^^ + >(r)^ = P, + P™. (5) 

Pj o 

where i?g = ks/rriu is the gas constant, is the atomic 
mass unit, fcs is the Boltzmann constant, and the mean 
atomic weight p. = (2X -I- 3/4Y)-^ = 0.588. We take 
the H mass fraction X = 0.76 and the He mass fraction 
Y = 0.24. In the resuhing models T >10,000 K except 
near the very surface, so the approximation for fi assumes 
that H and He are fully ionized. In the final models 
with masses near 800 M© the radiation pressure is of 
considerable importance. 

Given a mass M and an estimate for the outer ra- 
dius i?,, the hydrostatic structure is integrated outward 
from the center. At each point p{r) and T(r) are used 
to determine the Rosseland mean opacity k from a zero- 
metallicity table from OPAL (Iglesias & Rogers 1996), 



supplemented by a table from Lenzuni et al. (1991) for 
T < 6000 K. The photosphere is defined by the hydro- 
static condition 

kP = (6) 

where g is the acceleration of gravity. This point cor- 
responds to inward integrated optical depth r « 2/3. 
When it is reached, the local temperature is set to Toff 
and the stellar radiated luminosity is therefore 

= AnRlaBT^s (7) 

where i?* is the photospheric radius. 

The thermal equilibrium condition is L* = Ltot, where 
the total energy supply for the star (see below) is domi- 
nated during the DS phase by the DM luminosity Ljjm- 
Thus the next step is to determine L_dm from equation 
([2]) and from the density distribution of dark matter (dis- 
cussed in the next subsection). Depending on the value 
of LoM compared to L*, the radius is adjusted, and the 
polytrope is recalculated. The dark matter heating is 
also recalculated based on the revised baryon density dis- 
tribution. The radius is iterated until the condition of 
thermal equilibrium is met. 

2.2. DM Densities 

The value of the DM density inside the star is a key 
ingredient in DM heating. We start with a 10® Mq halo 
composed of 85% DM and 15% baryons. We take an ini- 
tial Navarro, Frenk, & White profile (1996; NFW) with 
a concentration parameter c = 2 at z = 20 in a standard 
ACDM universe. We note at the outset that our results 
are not dependent on high central DM densities of NFW 
haloes. Even if we take the initial central density to be a 
core (not rising towards the center of the halo), we obtain 
qualitatively the same results of this paper; we showed 
this in Freese et al (2008b). 

Enhanced DM density due to adiabatic contraction: 
Further density enhancement is required inside the star 
in order for annihilation to play any role. Paper I rec- 
ognized a key effect that increases the DM density: adi- 
abatic contraction (AC). As the gas falls into the star, 
the DM is gravitationally pulled along with it. Given the 
initial NFW profile, we follow its response to the chang- 
ing baryonic gravitational potential as the gas condenses. 
Paper I used a simple Blumenthal method, which as- 
sumes circular particle orbits (Blumenthal et al. 1986; 
Barnes & White 1984; Ryden & Gunn 1987) to ob- 
tain estimates of the density. Subsequently Freese et 
al. (2008b) did an exact calculation using the Young 
method (Young 1980) which includes radial orbits, and 
confirmed our original results (within a factor of two). 
Thus we feel confident that we may use the simple Blu- 
menthal method in our work. We found 

p^~5(GeV/cm3)(n,/cm3)0-«i, (8) 

where is the gas density. For example, due to this 
contraction, at a hydrogen density of 10^'^/cm'^, the DM 
density is 10^^ GeV/cm"^. Without adiabatic contrac- 
tion, DM heating in the first stars would be so small as 
to be irrelevant. 

Enhanced DM due to accretion: As further gas accretes 
onto the DS, driving its mass up from ~ IMq to many 
hundreds of solar masses, more DM is pulled along with 
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it into the star. At each step in the accretion process, 
we compute the resuhant DM profile in the dark star by 
using the Blumenthal et al. (1986) prescription for adi- 
abatic contraction. The DM density profile is calculated 
at each iteration of the stellar structure, so that the DM 
luminosity can be determined. 

The continued accretion of DM sustains the dark star 
for ~ 10^ years. Then, unless the DM is supplemented 
further, it runs out, that is, it annihilates faster than it 
can be resupplied. At this point the IO^Mq DM halo has 
a tiny "hole" in the middle that is devoid of DM. One 
might hope to rcpopulatc the DM by refilling the "hole," 
but this requires further study. On a conservative note, 
this paper disregards this possibility. Without any fur- 
ther DM coming into the star (which is now very large), 
it would contract and heat up to the point where fusion 
sets in; then the star joins the Main Sequence (MS). This 
is one possible scenario which we examine in the paper. 
Alternatively, there is another possible mechanism to re- 
plenish the DM in the dark star that may keep it living 
longer: capture. 

Enhanced DM due to capf,ure: At the later stages, the 
DM inside the DS may be rcpopulated due to capture 
(locco 2008; Freese et al. 2008c). This mechanism re- 
lies on an additional piece of particle physics, scattering 
of DM particles off the nuclei in the star. Some of the 
DM particles bound to the IO^Mq DM halo occasionally 
pass through the DS, and if they lose enough energy by 
scattering, can be trapped in the star. Then more DM 
would reside inside the DS and could continue to provide 
a heat source. This capture process is irrelevant during 
the initial stages of the DS, and only becomes impor- 
tant once the gas density is high enough to provide a 
significant scattering rate, which happens when the DS 
is already quite large (> lOOM©). For capture to be sig- 
nificant, there are two assumptions: (1) the scattering 
cross section must be high enough and (2) the ambient 
density of DM must be high enough. Regarding the first 
assumption: Whereas the existence of DS relies only on 
a standard annihilation cross section (fixed by the DM 
relic density), the scattering cross section is more specu- 
lative. It can vary over many orders of magnitude, and 
is constrained only by experimental bounds from direct 
detection experiments (numbers discussed below). Re- 
garding the second assumption: The ambient density 
is unknown, but we can again obtain estimates using 
the adiabatic contraction method of Paper 1. Indeed it 
seems reasonable for the ambient DM density to be high 
for timescales up to at most tens-hundreds of millions of 
years, but after that the lO^M© host for the dark star 
surely has merged with other objects and the central re- 
gions are likely to be disrupted; simulations have not yet 
resolved this question. Then it becomes more difficult for 
the DS to remain in a high DM density feeding ground 
and it may run out of fuel. 

We therefore consider two different situations: adia- 
batic contraction with and without capture; i.e., with 
and without additional capture of DM at the later stages. 
In the latter case we focus on "minimal capture," where 
the stellar luminosity has equal contributions from DM 
heating and from fusion, as described further below. In 
both cases, we find the same result: the dark stars con- 
tinue to accrete until they are very heavy. 



2.3. Evolution 

The initial condition for our simulations is a DS of 3 

Mq, roughly the mass where the assumption of nearly 
complete ionization first holds. We find an equilibrium 
solution for this mass, then build up the mass of the star 
in increments of one solar mass. 

We allow surrounding matter from the original bary- 
onic core to accrete onto the DS, with two different as- 
sumptions for the mass accretion: first, the variable rate 
from Tan & McKce (2004) and second, the variable rate 
from O'Shea & Norman (2007). The Tan/McKee rate 
decreases from 1.5 x 10~^ M0/yr at a DS mass of 3 Mq 
to 1.5 X 10-3 Mg/yr at 1000 Mq. The O'Shea/Norman 
rate decreases from 3 x 10"^ Mq/jt at a DS mass of 3 
Mq to 3.3 X 10"'* Mo/yr at 10l}0 Mq. The accretion 
luminosity arising from infall of material from the pri- 
mordial core onto the DS is not included in the radiated 
energy. It is implicitly assumed to be radiated away in 
material outside the star, for example an accretion disk 
(McKee & Tan 2008), and in any case, at the large radii 
involved here, this energy is small compared with the 
DM annihilation energy. 

We remove the amount of DM that has annihilated at 
each stage at each radius. We continue stepping up in 
mass, in increments of 1 Mq, until we reach 1000 Mq, the 
Jeans mass of the core (Bromm & Larson 2004) or until 
all of the following conditions are met: 1) the total mass 
is less than 1000 Mq, 2) the dark matter available from 
adiabatic contraction is used up; 3) nuclear burning has 
set in, and 4) the energy generation from gravitational 
contraction is less than 1% of the total luminosity of the 
star. The last two of these conditions define the zero-age 
main sequence (ZAMS). The time step is adjusted for 
each model to give the correct assumed accretion rate. 
The final mass can be less than 1000 Mq because of 
feedback effects of the ionizing stellar radiation, which 
can reverse the infall of accreting gas (McKee & Tan 
2008). The feedback effect is included in an approximate 
manner, starting when the surface temperature reaches 
50,000 K. 

At the beginning of the evolution, the Schwarzschild 
stability criterion for convection shows that the model is 
fully convective; thus the polytrope of n = 1.5 is a good 
approximation to the structure. However as more mass is 
added and internal temperatures increase, opacities drop 
and, once the model has reached about 100 Mq, radia- 
tive zones start to appear at the inner radii and move 
outward. Beyond about 300 Mq the model is almost 
fully radiative. The evolutionary calculation takes into 
account this effect approximately, by gradually shifting 
the value of n up from 1.5 to 3 in the appropriate mass 
range. In the later stages of evolution, the constant value 
n = 3 is used. 

2.4. Energy Supply 

The energy supply for the star changes with time and 
comes from four major sources: 

The general thermal equilibrium condition is then 

L*=itot. (10) 

The contribution Ldm from adiabatically contracted 
DM dominates the heat supply first. Its contribution 
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been discussed previously and is given in equation (p|). 
Later, once the adiabatically contracted DM becomes 
sparse, the other effects kick in: gravity, nuclear burning, 
and captured DM. We discuss each of these contributions 
in turn. 

2.4.1. Gravity 

At the later stages of evolution, when the supply of 
DM starts to run out, the star adjusts to the reduction 
in heat supply by contracting. For a polytrope of n = 3, 
the total energy of the star is 

3 GAP 

Etot = —7: ^ Eth + ^-rad (11) 

z H 

where Eth is the total thermal energy and the radiation 
energy i?rad — J aT'^dV . Using the virial theorem 



3GA/2 
2 R 



2E< 



th 



E, 



ad 



0, 



(12) 



we eliminate i^th from equation (jlip and obtain the grav- 
itational contribution to the radiated luminosity 
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2.4.2. Nuclear Fusion 



Once the star has contracted to sufficiently high in- 
ternal temperatures ^ lO^K, nuclear burning sets in. 
For a metal-free star, we include the following three en- 
ergy sources: 1) the burning of the primordial deuterium 
(mass fraction 4x 10~^) at temperatures around 10^ K, 2) 
the equilibrium proton-proton cycle for hydrogen burn- 
ing, and 3) the triple-alpha reaction for helium burning. 
At the high internal temperatures (2.75 x 10^ K) needed 
for the star to reach the ZAMS, helium burning is an im- 
portant contributor to the energy generation, along with 
the proton-proton cycle. We stop the evolution of the 
star when it reaches the ZAMS, so changes in abundance 
of H and He are not considered. Then Lnuc — J ^nucdAI 
where enuc is the energy generation rate in erg g~^ s^^; it 
is obtained by the standard methods described in Clay- 
ton (1968). The astrophysical cross section factors for the 
proton-proton reactions are taken from Bahcall (1989, 
Table 3.2), and the helium-burning parameters are taken 
from Kippenhahn & Weigert (1990). 

2.4.3. Feedback 

During the initial annihilation phase, the surface tem- 
perature Teff remains at or below 10^ K. However once 
gravitational contraction sets in, it increases substan- 
tially, reaching « 10^ K when nuclear burning starts. 
In this temperature range, feedback effects from the stel- 
lar radiation acting on the infalling material can shut off 
accretion (McKee & Tan 2008). These effects include 
photodissociation of H2, Lyman a radiation pressure, 
formation and expansion of an HH region, and photo- 
evaporation of a disk. These effects become important 
at a mass of « 50 M0 in standard Pop HI calculations of 
accretion and evolution, but the effect is suppressed un- 
til much higher masses in the DM case because of much 
lower Toff. Noting that 50 M0 on the standard metal-free 
ZAMS corresponds to T^s « 50,000 K (Schaerer 2002), 
we apply a linear reduction factor to the accretion rate 



above that temperature, such that accretion is shut off 
completely when T^s =100,000 K. This cutoff generally 
determines the final mass of the star when it reaches the 
ZAMS. 

2.4.4. Capture 

The contribution icap is determined by calculating the 
annihilation rate of captured DM 

(14) 
(15) 



Eca,p — 2772T^rcap, 



where 



d xp {(jv{r))/m^ 



and Pca.pi'r) is the density profile of captured DM. DM 
heating from captured DM has previously been discussed 
in detail by Freese et al. (2008c) and locco (2008). 

The density profile of the captured DM has quite a dif- 
ferent shape than the adiabatically contracted (AC) DM 
density profile discussed previously. Whereas the AC 
profile is more or less constant throughout most of the 
star (see Fig. 3), the captured DM profile is highly con- 
centrated toward the center with an exponential falloff 
in radius, as we will now show. Once DM is captured, 
it scatters multiple times, thermalizes with the star, and 
sinks to the center of the star. The DM then assumes the 
form of a thermal distribution in the gravitational well 
of the star 



Pcap 



/Eq/ThPo,ce' 
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Here Tc, the star's central temperature, is used since the 
DM has the same temperature as the stellar core; 0(r) is 
the star's gravitational potential; po,c is the central den- 
sity of the captured DM; and /xh and /eq are discussed 
below. One can immediately see another difference with 
the AC profile: for captured DM the profile is more cen- 
trally concentrated for larger DM particle mass, whereas 
the density profile for adiabatically contracted DM is in- 
dependent of particle mass. 

The quantity /xh is the probability that a captured 
DM particle has thermalized with the star. 



/Th = 1 



5-ti/TTh 



(17) 



ti is the time since a DM particle has been captured and 
Txh is the thermalization time scale of captured DM. 

Captured DM does not immediately thcrmalizc with 
the star; it requires multiple scatters. The thermalization 
time scale for the DM is 



2ascV 



(18) 



as previously discussed in Freese et al. (2GG8c) . Here Vcsc 
is the escape velocity from the surface of the star, uh is 
the average stellar density of hydrogen, and toh is the 
proton mass. We have only considered spin-dependent 
scattering, which includes scattering off hydrogen with 
spin 1/2, but not off Helium which has spin 0. 

Heating from captured DM is unimportant until the 
star approaches the MS; the thermalization time scale 
is very long (millions of years) and /xh ~ 0. Once the 
DM supplied from adiabatic contraction runs out, the 
star begins to contract, and the thermalization time scale 
becomes very short, on the order of a year. Thus, the star 
rapidly transitions from /xh ~ to /xh ~ 1. The DM is 
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either in thermal equihbrium or it is not. With this in 
mind, ti has been set to 100 yrs for all of the simulations. 
Only TTh is recalculated at every time step, which is a 
reasonable prescription; ti set between a few thousand 
years and a few years makes no quantitative difference 
upon the evolution of the DS. 

The ratio /eq gives the ratio of DM particles currently 
inside of the star due to capture, to the maximum 
number of DM particles that could possibly be present. 
No, 

/E.^£. (19) 

At the start of the simulation, the DS is too diffuse to 
capture any DM and /sq ~ 0. Once the star goes onto 
the main sequence, the maximum number is reached and 

/Eq « 1. 

We can now give an explicit expression of f^q by solv- 
ing for the total number of particles N bound to the 
star. This number is determined by a competition be- 
tween capture and annihilation, 

iV = C-2r, (20) 

where 'dot' refers to differentiation in time. The total 
number of DM particles captured per unit time C was 
calculated following Gould (1987a). The total number 
of DM particles annihilated per unit time F contains a 
factor of two since two particles are annihilated per an- 
nihilation.). Equation (|^D|) can be solved for N, 

N = Not&n\i{t2/Teq) (21) 

where No = \/ C/Ca, and r^q — I/^/CUa, and Ca = 
2r/iV2 (Griest & Seckel 1987). 

When the annihilation rate equals the capture rate, N 
goes to No- This transition occurs on the equilibrium 
time scale Teq. The equilibrium time scale is recalculated 
at each time step. The time t2 in equation (|21|) counts 
the time since capture is turned on. In equation (|16p . we 
now set 

N 

/Eq = -r^ - i&iMh/Teq). (22) 

Conceptually ^2 is not the same as ti but at a practical 
level it is reasonable to set t2 — ti. In the first place, ti 
refers to the time since a single DM particle has been 
captured and in fact ti is different for each captured DM 
particle. On the other hand, ^2 is a single time for all 
captured and thermalizcd DM and starts when capture 
turns on. Regardless, once the star descends towards the 
main sequence, /Eq and /xh both go to unity. The con- 
nection is not a coincidence. First r^q is a factor of a 
few to a hundred less than TTh throughout the simula- 
tion. In addition, t2 > ti except for the very first DM 
particle captured. Fixing t2 = ii does not change the 
behavior of /Eq in relation /xh as just described. Hence, 
for simplicity the simulation assumes ti = t2. 

In the limit that /sq and /xh go to unity, then Tcap = 
C/2. Hence, Leap simplifies to 

2 

where the 2/3 accounts for the loss of neutrinos as pre- 
viously discussed in connection with equation Also 
roughly Leap oc p^asc and is independent of m^. Hence 



Leap scales with both the ambient background DM den- 
sity Px , and the scattering cross section of the DM off of 
baryons CTsc- 

Physically reasonable values for the scattering cross 
section and background DM densities imply that Leap 
can affect the first stars once they approach the MS. Di- 
rect and indirect detection experiments require Usc < 
4 X 10-39 (Savage et al. 2004; Chang et al. 2008; Sav- 
age et al 2008; Hooper et al. 2009); a WIMP's a^c can 
be many orders of magnitude smaller than the exper- 
imental bound. The ambient background DM density, 
p^, can be very high, on the order of 10^'*(GeV/cm3) 
(as obtained from Eq. (8)) due to adiabatic contraction. 
To demonstrate the "minimal" effect of capture, we set 
Px = 1.42 X 10i''(GeV/cm3) and Usc = lO'^^cm^, chosen 
so that about half of the luminosity is from DM capture 
and the other half from fusion for the standard 100 GeV 
case at the ZAMS. These parameters have been used in 
all cases where capture has been turned on. Even min- 
imal capture shows that DM capture lowers the star's 
surface temperature and increases the mass of the first 
stars. 

We note that the product of ambient density and cross 
section could be very different than the values used in 
this paper. If this product is much smaller than the val- 
ues considered in the 'minimal capture' case, then fusion 
completely dominates over capture and one may ignore 
it. On the other hand, if this product is much larger 
than the 'minimal capture' case, then fusion is negligible 
and the DS may continue to grow by accretion and be- 
come much larger than discussed in this paper; this case 
will be discussed elsewhere. The 'minimal' case is thus 
the borderline case between these two possibilities, and 
lasts only as long as the DS remains embedded in a high 
ambient DM density assumed here. 

3. RESULTS 

Here we present our results for the evolution of a 
DS, starting from a 3Mq star powered by DM annihi- 
lation (defined as t = 0), building up the stellar mass 
by accretion, and watching as the other possible heat 
sources become important. We have considered a range 
of DM masses: 1 GeV, 10 GeV, 100 GeV, 1 TeV and 10 
TeV, which characterizes the plausible range for WIMP 
masses. Since the DM heating scales as {av)/mx, this 
range of particle masses also corresponds to a range of 
annihilation cross sections. To account for the uncer- 
tainties associated with the accretion rate, the study has 
used two different accretion rates: one from Tan & Mc- 
Kee (2004) and the other from O'Shea & Norman (2007) 
(as discussed in section 12.31 above). We consider the 
case of (i) no capture and (ii) "minimal capture" (de- 
fined as above), with p^ = 1.42 x 10^'^(GeV/cm3) and 
(Jsc = lO-'^^cm^, chosen so that about half of the lumi- 
nosity is from DM capture and the other half from fusion 
for the standard 100 GeV case at the ZAMS. In this pa- 
per we will not treat the case of maximal capture, where 
the DM heating would dominate (over all other sources 
including fusion) so that the DS would continue to grow 
as long as it is fed more DM . 

We here show the results of our simulations. First we 
briefly outline our results and present our tables and fig- 
ures. These are then discussed in detail, case by case, in 
the remainder of the results section. 
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TABLE 1 

Properties and Evolution of Dark Stars for = 1 GeV, M from Tan & McKee, (av) = 3 x 10~^® cm^/s. For the lower 

STELLAR masses, IDENTICAL RESULTS ARE OBTAINED FOR THE CASES OF NO CAPTURE AND MINIMAL CAPTURE; HENCE THE FIRST FOUR ROWS 
APPLY EQUALLY TO BOTH CASES. THE FIFTH AND SIXTH ROWS APPLY ONLY TO THE NO CAPTURE CASE WHILE THE LAST TWO ROWS APPLY 

ONLY WHEN CAPTURE IS INCLUDED. 



M, 


L, 


R, 




Pc 


Tc 


Mdm 


Po.DM 


t 


(Me) 


(lO'Lo) 


(lO^^cm) 


(lO^K) 


(gm/cm^) 


(lO'^K) 


(lO^^g) 


(gm/cm'') 


(lO'yr) 


106 


9.0 


240 


5.4 


1.8 X 10"* 


0.1 


36 


4.5 X 10"^^ 


19 


371 


17 


270 


5.9 


5.2 X 10"*^ 


0.20 


64 


2.0 X 10"^^ 


119 


690 


6.0 


110 


7.5 


5.9 X 10"^ 


1.0 


10 


6.8 X 10"" 


280 


756 


3.3 


37 


10 


2.2 X 10"" 


3.4 


1.2 


4.8 X 10"^° 


310 


793 


4.5 


5.1 


31 


8.9 X 10"^ 


25 


^ 0.0 


^ 0.0 


330 


820 


8.4 


0.55 


111 


115 


276 


« 0.0 


« 0.0 


430 


Alternative properties 


















for final two masses 


















with capture included: 


















793 


4.6 


5.7 


30 


7.1 X 10"^ 


24 


8.4 X 10"" 


2.94 X 10"** 


328 


824 


8.8 


.58 


109 


102 


265 


8.1 X 10"^ 


2.14 X 10"** 


459 



TABLE 2 

Properties and Evolution of Dark Stars for = 100 GeV, M from Tan & McKee, (av) = 3 x 10"^^ cm^/s. Entries as 

DESCRIBED IN TABLE 1. 



(Mo) 


L, 

(lO^Lo) 


R, 

(lO^^cm) 


Tcs 
(lO^K) 


pc 

(gm/cm^) 


(lO^K) 


Mdm 
(lO^ig) 


Po,DM 

(gm/cm*) 


t 

(lO^r) 


106 


1.0 


70 


5.8 


1.1 X 10"** 


0.4 


16 


1.8 X 10"'' 


19 


470 


4.8 


84 


7.8 


2.0 X 10"" 


1.4 


41 


2.0 X 10"'' 


171 


GOO 


5.0 


71 


8.5 


i.l X 10 ' 


1.8 


:!() 


1.9 X 10 '' 


2:!5 


716 


6.7 


11 


23 


2.0 X lO""" 


15 


2.7 


3.8 X lO"*" 


303 


756 


6.9 


2.0 


56 


2.8 


78 


^ 0.0 


^ 0.0 


330 


779 


8.i 


().."•)."•) 


ill) 


121) 


281) 


=3 0.0 


=s 0.0 


:i87 


Alternative properties 
for final two masses 
with capture included: 




756 (c) 


6.8 


2.0 


55 


2.55 


76 


7.5 X lO"-""' 


5.3 X lO"-""' 


328 


787(c) 


8.5 


0.58 


108 


105 


270 


2.2 X 10"^ 


5.3 X 10"" 


414 



TABLE 3 

Properties and Evolution of Dark Stars for = 10 TeV, M from Tan & McKee, (crv) = 3 x 10"^® cm^/s. Entries as 

DESCRIBED IN TABLE 1. 





L, 


R, 




Pc 


Tc 


Mdm 


Po.DM 


t 


(Mq) 


(10«Lq) 


(lO^^cm) 


(lO^K) 


(gm/cm^) 


(lO^K) 


(lO^ig) 


(gm/cm^) 


(lo'yr) 


106 


0.1 


22 


6.0 


3.8 X 10"^ 


1.3 


7.8 


1.8 X 10"* 


19 


310 


0.59 


21 


9.3 


4.4 X 10"" 


3.6 


15 


1.5 X 10"^ 


88 


399 


3.3 


6.6 


25 


3.6 X 10"^ 


16 


6.6 


9.9 X lO"'^ 


134 


479 


3.0 


2.9 


32 


0.5 


40 


2.1 


1.1 X 10"** 


172 


552 


5.2 


0.56 


97 


79 


220 


^ 0.0 


fa 0.0 


230 


552 


5.3 


0.47 


107 


136 


270 


« 0.0 


« 0.0 


256 


Alternative properties 


















for final two masses 


















with capture included: 


















550 (c) 


5.2 


0.6 


95 


66 


210 


8.5 X 10"'^ 


5.3 X 10" 


230 


553 (c) 


5.4 


0.48 


105 


127 


260 


5.7 X 10"'^ 


1.2 X 10^ 


266 
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Tables 1-3, each for a fixed m^, present the evolu- 
tion of some basic parameters of the star throughout the 
DS phase: stellar mass M*, luminosity L*, photospheric 
radius i?,, surface temperature Tcs, the star's baryonic 
central density pc, central temperature Tc, the amount 
of DM inside of the star A/dm, the DM's central density 
Po,DM, as well as the elapsed time t since the start of the 
calculation. Mdm and Po,DM give the mass and density 
of DM either from adiabatic contraction (the first 6 rows) 
or from capture (the last two rows). 

In the tables, the stellar masses were chosen to roughly 
characterize important transitions for the DS. The first 
mass was chosen at roughly the time when a dark star 
begins to take on the size, temperature, and luminos- 
ity which are characteristic of the DS phase, which is 
around 100 Mq. The second mass is given when the DS 
acquires its maximum amount of DM, which varies from 
case to case. The third mass is towards the end of the DS 
adiabatic phase. The fourth row characterizes the star's 
properties just prior to running out of adiabatically con- 
tracted DM. The fifth row is given just after the star runs 
out of adiabatically contracted DM, while the sixth gives 
the star's properties on the ZAMS. Capture only becomes 
important once the adiabatically contracted DM has run 
out. Prior to running out, the stellar evolution is the 
same with or without DM capture. For the case with 
capture included, the seventh and eighth row should be 
used instead of the fifth and sixth row. 

Further features of the calculations are shown in Fig- 
ures 1-4, which give, respectively, (1) the evolution of 
all cases without capture in the Hertzsprung-RusscU dia- 
gram, (2) the luminosity as a function of time for the 100 
GeV case with and without capture, (3a) the evolution of 
the density distribution of the baryonic and dark matter 
for the 100 GeV case without capture, (3b) the distribu- 
tion of dark matter heating rate for the 100 GeV case 
with capture, and (4) the evolution of various quantities 
for all cases without capture. 

3.1. Canonical Case without Capture 

Our canonical case takes = 100 GeV, the accretion 
rate as given by Tan & McKee (2004), and {av) = 3 x 
10~^^ cm'^s"^ for the calculation of the annihilation rate. 

The initial model (defined to be t = 0) has a mass of 
3 Mq with a surface temperature of 4100 K, a radius of 
3 X 10^'^ cm, a luminosity of 4.5 x 10"* Lq and a fully 
convective structure. The DS stays fully convective for 
stellar masses below 50 M© at which point a small ra- 
diative zone appears. The star makes a transition from 
convective to radiative in the M, = (100 — 200) Mq mass 
range, with a radiative zone growing outward from the 
center, and then becomes fully radiative (but for a small 
convective region at the surface) for M, > 300 M© . The 
star stays almost fully radiative for the remainder of the 
DS phase. Above 750 M©, hydrogen and helium burning 
set in, and a new convective zone appears at the center 
of the star. 

Evolution of Dark Matter Density: The DM densities 
inside of the star change due to annihilation and the evo- 
lution of the dark star's baryonic structure. Figure [3] 
plots the baryonic and DM density profiles for the stan- 
dard case at several different stellar masses. As the star's 
mass grows, the DM and baryonic densities also increase. 
The DM profile for the last model, which has a stellar 



mass of 779 M©, is not plotted because at this point all 
of the adiabatically contracted DM has been used up. 

The adiabatically contracted DM density profile is 
more or less constant with radius and then falls off to- 
ward the surface of the star. The DM densities for the 
earlier models show a density increase toward the cen- 
ter of the star due to the original NFW profile, e.g. in 
the 300 M© model. The slight peak flattens out at higher 
stellar masses as can be seen in the 500 M© models, since 
regions with a higher DM density annihilate more rapidly 
than regions with a lower DM density. Once the peak is 
worn away, the central DM profile stays flat and ulti- 
mately decreases. 

The DM and baryonic densities evolve slowly for much 
of the DS phase. In fact between 100 M© and 600 M©, 
the DM central density barely changes (Table 2). Ini- 
tially, the DM mass grows (Fig. 4e) and reaches a max- 
imum when the total mass is about 450 M©, at which 
point the star has about 0.2 M© of DM. At this point, 
the DM heating becomes insufficient to support the star. 
The star begins to slowly contract and boosts the DM 
density, but the boost is still modest up to a time of 
2 X 10^ years or until the star reaches 600 M©. 

Aroimd 600 M©, the star begins to rapidly shrink. 
Note the rapid dominance of gravitational energy at this 
time (Fig. [2]). The DM densities shoot up, but the total 
amount of DM inside of the star precipitously drops. By 
750 M© , the DM has run out and the star begins to tran- 
sition into a standard Pop. Ill star. The star reaches the 
MS at around 780 M© at which point feedback effects 
cut off any further accretion. 

Evolution of temperature and energy generation: The 
surface properties of the DS evolve rapidly between the 
initial model at 3 M© and 100 M©. The luminosity 
shoots up by a factor of 25, the radius by a factor of 
3, and the star's surface temperature goes from 4100 K 
to 5800 K. 

The star looks radically different from a metal-free star 
on the MS. Even at 100 M© , the DS is stih very cool with 
a surface temperature of around 5800 K vs 100, 000 K (for 
the standard ZAMS) and has a central temperature of 
only about a half a million degrees Kelvin vs 200 million 
degrees (ZAMS). The star is also much more extended, 
« 7 X 10^3 cm vs « 5 X 10" cm (ZAMS). 

Between 100 — 600 M©, the star's luminosity evolves 
less rapidly, increasing by a factor of 5 to a value of 
5 X 10^ L©. The star's surface temperature reaches 
7800 K, while the radius remains almost constant, pass- 
ing through a broad maximum at about 8.5 x 10^'^ cm. 
The central temperature is still quite cool (1.8 x 10^ K) 
compared to a metal-free ZAMS star. Most of the star is 
still too cool to even burn deuterium and lithium. DM 
continues to power the star, as is apparent in Figure [2l 

Above 600 M©, the DS begins to transform into a stan- 
dard metal- free star. As the DM begins to run out the 
surface temperature begins to shoot up and the stellar ra- 
dius shrinks. Just prior to running out of DM, at 716 M©, 
the star's surface reaches 23, 000 K, the radius has shrunk 
to 10^^ cm and the central temperature has shot up to 15 
million K. Near the beginning of this rapid contraction 
the star is sufficiently hot to burn deuterium rapidly, as 
can be seen from the spike of nuclear heating in Figure [2] 
just before 0.3 Myr; this spike lasts about 50,000 years. 
The burning of deuterium, however, remains an insignifi- 
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Fig. 1. — Evolution in the H-R diagram, for 5 different particle masses as indicated in the legend. The luminosity is given in solar units 
and the temperature is in Kelvin. The dots indicate a series of time points, which are the same for all cases. The open dots distinguish the 
100 GeV (B) case from the 100 GeV case. All cases are calculated with the accretion rate given by Tan & McKee (2004) except the curve 
labelled 100 GeV (B), which is calculated with the rate given by O'Shea & Norman (2007). The metal-free zero-age main sequence (Z = 
ZAMS) is taken from Schaerer (2002). The peak below lO'^K in the 1 GeV case is due to the overwhelming DM luminosity in this case. 



cant power source for the star. As the DS runs out of DM 
gravity and nuclear burning begin to turn on as energy 
sources for the star. Just at the point where all three en- 
ergy sources (DM heating, gravity, and nuclear burning) 
all contribute simultaneously, the luminosity increases by 
about 50 % and then decreases again; one could call this 
a "flash" which lasts about 20,000 years and is the "last 
gasp" of adiabatically contracted DM annihilation. 

Near 750 Mq the energy generation in the star becomes 
dominated by, first, gravity, and then nuclear burning. 
The star's surface temperature is just above 50,000 K 
(Table 2) and feedback effects have started to cut off the 
accretion of baryons. Eventually, the star begins to settle 
onto the MS and nuclear burning becomes the dominant 
power source. The surface temperature has shot up to 
109,000 K, and finally feedback effects have stopped any 
further accretion. On the MS, the radius is 5.5 x 10^^ cm 
and the central temperature has gone up to 280 million 
K. At the final cutoff, which is our zero- age main se- 
quence model, the star has a mass of 779 M© . 

The final model on the metal-free MS compares well 
with that produced by a full stellar structure code 
(Schaerer 2002). Schaerer 's radius is slightly larger at 
8.4 X 10^^ cm vs. 5.5 x 10^-^ cm for the polytrope. The 



polytrope's surface temperature is slightly too high at 
109,000 K vs 106,000 K and the polytrope's luminosity 
is too low by a factor of about 2. The difference is not 
surprising; we have not done the radiative energy trans- 
port. Remarkably, the numbers are quite close despite 
the polytropic assumption. 

3.2. General Case without Capture 

In addition to the canonical case, we have studied many 
other parameter choices as described previously: a vari- 
ety of WIMP masses (or, equivalently, annihilation cross 
sections) and a variety of accretion rates. In this sec- 
tion we still restrict ourselves to the case of no capture. 
In all cases, we obtain the same essential results. While 
one can see in Figure [3] that there is some scatter of the 
basic properties of the star, yet we find the following 
general results for all cases: DS evolution occurs over an 
extended period of time ranging from (2 — 5) x 10^ years. 
All cases have an extended radius (2 x 10^^ — 3 x 10^^)cm, 
are cool {Tcb < 10,000 K) for much of the evolution, and 
have a large final mass (500 - 1000) Mq. All of the DS 
cases studied are physically very distinct from standard 
Pop. Ill stars. The standard Pop. Ill star is physically 
more compact with a radius of only w 5 x lO^^cm, and 
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Fig. 2. — Luminosity evolution for the 100 GeV case as a function 
of time {lower scale) and stellar mass {upper scale). The solid 
(red) top curve is the total luminosity. The lower curves give the 
partial contributions of different sources of energy powering the 
star a) {upper frame) without capture, and b) {lower frame) with 
'minimal' capture. In both frames, the total luminosity is initially 
dominated by DM annihilation (the total and annihilation curves 
are indistinguishable until about 0.3 Myr after the beginning of the 
simulation); then gravity dominates, followed by nuclear fusion. In 
the lower frame, capture becomes important at late times. 

hotter with a surface temperature of «100,000 K. The 
estimated final mass for the standard Pop III.l star is 
also smaller, w 100 Mq, because the high surface tem- 
perature leads to radiative feedback effects that prevent 
accretion beyond that point (McKee & Tan 2008). 

Figure[T]shows the evolution curves in the H-R diagram 
for a variety of cases. The figure shows the evolution for 
various particle masses with the Tan/McKee accretion 
rate, and compares the 100 GeV case with a calculation 
using the O 'Shea/Norman accretion rate. 

The DS evolution curves on the HR diagram are qual- 
itatively similar for all of the different cases. The DS 
spends most of its lifetime at low temperatures below 
10,000 K. Eventually, DM begins to run out, the star 
contracts and heats up, and gravity and nuclear burn- 
ing begin to turn on. Just at the point where all three 
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Fig. 3. — a) (upper frame): Evolution of a dark star for m^,,. = 100 
GeV as mass is accreted onto the initial protostellar core of 3 Mq 
(for the case of no capture). The set of upper (lower) solid curves 
correspond to the baryonic (DM) density profile (values given on 
left axis) at different stellar masses and times. h){Lower frame): 
differential luminosity rdLYiyi/dr as a function of r for the masses 
indicated. 



power sources are important, the luminosity "flashes" 
(as seen in Fig. 2), i.e., increases by a 50 % for about 
20,000 years and then goes back down. Then, DM runs 
out. For the lower particle masses, as seen in Fig. 1, the 
luminosity drops slightly as gravity becomes the main 
and dominant power source. As the surface temperature 
continues to increase, nuclear burning becomes impor- 
tant when Teff ~ 30, 000 K. Above this temperature (see 
Fig. 1), the stellar luminosity again increases as the star 
approaches the ZAMS. Above 50,000 K, feedback effects 
turn on and begin to stop the accretion of baryons. Once 
the star reaches 100,000 K, accretion completely stops, 
and the star goes onto the ZAMS. 

While the basic DS evolution is similar for all cases, 
we now discuss the variation due to different parame- 
ters. The evolution curves differ depending upon the 
DM particle mass. For larger m^, the stellar mass and 
luminosity decrease such that the evolution curve on the 
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Fig. 4. — Dark Star evolution (for the case of no capture) illustrating dependence on DM particle mass and accretion rate(0'Shea/Norman 
vs. Tan/McKee). Different curves are the same as explained in the legend of Fig. 1. The quantities plotted are a) stellar mass (upper 
left), b) total radius (upper right), c) surface temperature (middle left), d) central stellar (baryon) density (middle right), e) amount of DM 
inside of the star (lower left), and f) radiative gradient at the center of the star (lower right). In the latter plot, central regions above the 
horizontal line are unstable to convection; hence this plot illustrates the onset of a central radiative zone. 



H-R diagram for the 1 TeV case is beneath that for the 
100 GeV case. The reduced luminosity is a consequence 
of the reduced energy production at the higher particle 
masses, since DM heating is inversely proportional to 
(see eq. ([T])). Then at higher particle masses a smaller 
radius is needed to balance the dark-matter energy pro- 
duction and the radiated luminosity. 

In the Tables and in Figure 3, one can see the den- 
sity profiles of the various components in the DS. The 
DM densities also depend upon m-^. As shown in the 
Tables, the average DM density in the star is an increas- 
ing function of m^. For example, at the time the DS 
reaches 100 Mq, the central DM densities vary between 
2 X 10~^ g/cm^ for = 10 TeV and 5 x 10"" g/cm^ 
for = 1 GeV. These densities are many orders of 
magnitude lower than the baryon densities in the DS, as 
can be seen in Figure 3. Here one can see "the power 
of darkness:" Despite the low DM densities in the star, 
DM annihilation is an incredibly efficient energy source 
since nearly all of the particle mass is converted to heat 
for the star. We also wish to point out that both baryon 
and DM densities in the DS arc far lower than in metal 
free main sequence stars of similar stellar mass; ZAMS 
stars have densities > 1 gm/cm'^. 

One can also see the dependence of the DS evolution 



upon accretion rate. Our canonical case, obtained us- 
ing the accretion rate from Tan & McKee (2004), can be 
contrasted with that using the O'Shea/Norman accretion 
rate. For 100 GeV particles, the DS evolution is identical 
for the two accretion rates until Toff ~ 7, 000 K (10^ yr). 
After that the luminosity of the O'Shea/Norman case 
shoots higher as a result of the higher accretion rate. 
Consequently the DS with the O'Shea/Norman accre- 
tion rate has a shorter lifetime than the DS with the 
Tan/McKee accretion rate (350 vs. 400 thousand years); 
a larger final mass (916Mo vs. TSOM©) and higher lu- 
minosity. 

The 1 GeV case displays some interesting behavior: 
one can see a spike in the evolutionary track on the H-R 
diagram. The DS reaches a maximum luminosity of 2 x 
10*" Lq and then subsequently drops by a factor of nearly 
ten, crossing beneath the 100 GeV case. Subsequently 
the 1 GeV curve then again crosses the 100 GeV case 
before settling onto the MS. The large early luminosity 
is due to extreme effectiveness of DM heating for the 1 
GeV case; DM heating is inversely proportional to to^ 
(eq. m- 

The stellar structure of all of the cases studied is similar 
to our canonical case, but the transition from a convec- 
tive to a radiative structure happens at different times, 
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depending upon the DM particle mass. The lower right 
panel in Fig. 4 plots the temperature gradient in the 
innermost few zones at the center of the DS, to show 
the onset of radiative energy transfer at the center of the 
star. Models above the line have convective cores while 
models below have radiative cores. As we increase the 
DM mass, the onset of radiative transport happens ear- 
lier, at lower stellar masses. For our canonical 100 GeV 
case, a small radiative zone appears when the mass of 
the DS reaches 50Mq. The = 1 GeV case, however, 
does not begin to form a radiative zone until the mass is 
nearly 400 Mq. The larger radii and cooler internal tem- 
peratures, compared with the canonical 100 GeV case, 
result in higher interior opacities, favoring convection. 

The timing of the transition from convective to radia- 
tive can be understood from the stellar structure. For 
greater than the canonical case, the DM heating is less 
effective since it is inversely proportional to m^. The 
star compensates by becoming more compact, thereby 
increasing the DM density and the internal temperature. 
In the relevant temperature range, the opacity scales as 
pT~^-^, and p scales as T^. The smaller radius actually 
results in a decrease in luminosity as well, so the radiative 
gradient 

dlogr ^ 3kLP 

d logP WnGacT'^Mr ^ ' 

is reduced, tending to suppress convection relatively 
early. Conversely if is less than lOOGeV, the star 
must be more puffy, which increases the radiative gradi- 
ent. Thus the transition will occur at a later stage. 

The cases studied have some differences in terms of 
lifetime, luminosity, and final mass. At the time the 
DS has grown to 100 M©, the 1 GeV case is an order 
of magnitude brighter than the 100 GeV case. Similarly, 
the 100 GeV case is about an order of magnitude brighter 
than the 10 TeV case. The DS lifetime grows shorter for 
increasing m^. The lifetime for the 1 GeV case is 4.3 x 10^ 
years while the lifetime for the 10 TeV case is only 2.6 x 
10^ years. The reduced amount of fuel available in the 10 
TeV case more than compensates for the lower luminosity 
in determining the lifetime. Also, the 1 GeV case has 
a higher final mass than the 10 TeV case (820 M© vs 
552 Mq); clearly a longer lifetime allows more mass to 
accrete. 

3.3. DM Evolution- Adiabatic Contraction Only 

Here we discuss the evolution of the DS once its DM 
fuel begins to run out. In this section, we do not consider 
capture; in the next section we will add the effects of 
capture as a mechanism to refuel the DM inside the star. 

Initially, as baryons accrete onto the DS, adiabatic con- 
traction pulls DM into the star at a sufficiently high rate 
to support the star, but eventually, adiabatic contraction 
fails to pull in enough DM at high enough densities to 
support it. 

Adiabatic contraction eventually fails for two main rea- 
sons. DM must be brought inside the star, first, to re- 
place the DM which has annihilated away and, second, 
to support an ever-increasing stellar luminosity as the 
stellar mass increases due to accretion. In fact, an ac- 
cretion rate that decreases as a function of time (Tan & 
McKee 2004) only compounds the problem, leading to a 
relatively shorter life span since even less DM is brought 



into the star. With an insuffient amount of DM at high 
enough densities, the star must contract to increase the 
DM heating rate to match the stellar luminosity. The 
contraction leads to higher surface temperatures such 
that feedback radiative effects eventually shut off accre- 
tion altogether. 

As the star contracts, a large fraction of DM previously 
inside of the star ends up outside of the star. The star's 
contraction amounts to at least two orders of magnitude 
in radius, leaving a large fraction of the DM outside the 
star. Thus DM annihilation accounts for only a fraction 
of the DM "lost" from the star. For instance in the 100 
GeV case, only a fifth of the DM is lost due to annihila- 
tion; the rest is "left outside". 

The "left outside" effect can be understood in a simple 
way. DM heating is mainly a volume effect as can be seen 
in Figure [21 Due to adiabatic contraction, poM does not 
scale perfectly with baryon density (eq. [8|). For instance, 
the total amount of DM inside of star of fixed mass (even 
without DM annihilation) decreases as the radius of the 
star shrinks. 

Without a way to replenish the DM such as with cap- 
ture, the DS will eventually run out of DM due to the 
combination of the "left outside" effect and DM annihi- 
lation. 

3.4. DS evolution with Capture included 

The DM inside the DS can be replenished by cap- 
ture. DM particles bound to the larger halo pass through 
the star, and some of them can be captured by los- 
ing energy in WIMP/nucleon scattering interactions. 
We remind the reader that we have defined our case 
of 'minimal capture' to have background DM density 
p^ — 1.42 X 10^°(GeV/cm'^) and scattering cross sec- 
tion ffsc = 10~'^^cm^ chosen so that about half of the 
luminosity is from DM capture and the other half from 
fusion for the standard 100 GeV case at the ZAMS. 

With capture, the stellar evolution is the same as in 
the cases previously discussed without capture through- 
out most of the history of the DS. The contribution to the 
heating from capture Leap only becomes important once 
DM from adiabatic contraction runs out. Once this hap- 
pens, gravity still initially dominates, but for the 'mini- 
mal capture' case we have considered, Leap is more im- 
portant than Lnuc until the star eventually reaches a new 
"dark" main sequence, where the two contributions be- 
come equal. During this phase the star is more extended 
and cooler than without capture, which allows for more 
baryons to accrete onto the star once the first dark star 
phase driven by adiabatic contraction ends. However for 
the standard minimal case the final mass with capture is 
at most 1% larger than that without. 

In a future paper we will consider the case of much 
higher capture rates than the minimal ones we have 
studied here. The capture rate depends on the ambi- 
ent DM density in which the DS is immersed and on the 
WIMP/nucleus scattering cross section. High ambient 
density could easily arise; in fact the adiabatic contrac- 
tion arguments applied to the DS from the beginning 
indicate that this is quite likely. On the other hand, the 
scattering cross section may be significantly lower than 
the current experimental bound adopted in our minimal 
capture case. In any case, with a sufficiently high capture 
rate, say two or more orders of magnitude higher than 
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the minimal value considered here, the star can stay DM 
powered and sufficiently cool such that baryons can in 
principle continue to accrete onto the star indefinitely, or 
at least until the star is disrupted. This latter case will 
be explored in a separate paper where it will be shown 
that the dark star could easily end up with a mass on the 
order of several tens of thousands of solar masses and a 
lifetime of least tens of millions of years. 

4. DISCUSSION 

We have followed the growth of equilibrium protostel- 
lar Dark Stars (DS), powered by Dark Matter (DM) an- 
nihilation, up to the point when the DS descends onto 
the main sequence. Nuclear burning, gravitational con- 
traction, and DM capture are also considered as energy 
sources. During the phase when dark-matter heating 
dominates, the objects have sizes of a few AU and cen- 
tral T « 10^ - 10^ K. Sufficient DM is brought into the 
star by contraction from the DM halo to result in a DS 
phase of at least several hundred thousand years (without 
DM capture). Because of the relatively low T^s (4000- 
10,000 K), feedback mechanisms for shutting off accre- 
tion of baryons, such as the formation of HII regions or 
the dissociation of infalling H2 by Lyman- Werner pho- 
tons, are not effective until DM begins to run out. The 
implication is that main-sequence stars of Pop. Ill are 
very massive. Regardless of uncertain parameters such 
as the DM particle mass, the accretion rate, and scatter- 
ing, DS are cool, massive, puffy and extended. The final 
masses lie in the range 500-1000 Mq, very weakly de- 
pendent on particle masses, which were assumed to vary 
over a factor of lO"*. 

One may ask how long the dark stars live. If there 
is no capture, they live until the DM they are able to 
pull in via adiabatic contraction runs out; the numerical 
results show lifetimes in the range 3 x 10^ to 5 x 10^ yr. 
If there is capture, they can continue to exist as long as 
they reside in a medium with a high enough density of 
dark matter to provide their entire energy by scattering, 
capture, and annihilation. 

In addition, the properties of the DM halo affect the 
nature of the DS. For instance, a larger concentration pa- 
rameter allows for more DM to be pulled in by adiabatic 
contraction. The additional DM prolongs the lifetime 
of the DS phase and increases the DS final mass. The 
canonical case with a concentration parameter c = 3.5 
rather than x = 2 has a final mass of 870 Mq vs. 779 Mq 
and has a hfetime of 442,000 yrs vs. 387,000 yrs. Fur- 
thermore, the concentration parameter of each DS halo 
will be different, depending upon the history of each halo. 
Thus the Initial Mass Function (IMF) of the ZAMS aris- 
ing from the growth of dark stars will be sensitive to the 
halo parameters. This effect will be studied thoroughly 
in a separate publication. 

DS shine with a few IO^Lq and one might hope to 
detect them, particularly those that live to the lowest 
redshift (z = in the most optimistic case). One may 
hope that the ones that form most recently are detectable 
by JWST or TMT and differentiable from the standard 
metal-free Pop. Ill objects. DS are also predicted to have 
atomic hydrogen lines originating in the warmer photo- 



spheres, and II2 lines arising from the infalling material, 
which is still relatively cool. 

The final fate of DS once the DM runs out is also dif- 
ferent from that of standard Pop III stars, again leading 
to different observable signatures. Standard Pop III stars 
are thought to be ~ 100 — 2OOAf0, whereas DS lead to 
far more massive MS stars. Heger & Woosley (2002) 
showed that for I4OM0 < M < 260Mq, pair instability 
supernovae lead to odd-even effects in the nuclei pro- 
duced, an effect that is not observed. Thus if Pop III.l 
stars are really in this mass range one would have to con- 
strain their abundance. Instead, we find far more mas- 
sive Dark Stars which would eventually become heavy 
main sequence (ZAMS) stars of 500-1000 Mq. Heger & 
Woosley find that these stars after relatively short life- 
times collapse to black holes. Or, if they rotate very 
rapidly, Ohkubo et al (2006) argue that they could be- 
come supernovae, leaving behind perhaps half their mass 
as black holes. In this case the presumed very bright su- 
pernova could possibly be observable. In either case, the 
elemental abundances arising from pair instability SN of 
standard Pop III (less massive) stars are avoided so that 
one might avoid the constraint on their numbers. Other 
constraints on DS will arise from cosmological consider- 
ations. A first study of their effects (and those of the 
resultant MS stars) on reionization have been done by 
Schleicher et al. (2008a, 2008b), and further work in this 
direction is warranted. 

The resultant black holes from DS could be very impor- 
tant. DS would make plausible precursors of the IO^Mq 
black holes observed at z = 6 (Li et al. 2007; Pelopessy 
et al. 2007); of Intermediate Mass Black Holes; of black 
holes at the centers of galaxies; and of the black holes 
recently inferred as an explanation of the extragalactic 
radio excess seen by the ARCADE experiment (Seiffert 
et al. 2009). However, see Alvarez et al. (2008) who 
present caveats regarding the growth of early black holes. 
In addition, the black hole remnants from DS could play 
a role in high-redshift gamma ray bursts thought to take 
place due to accretion onto early black holes (we thank 
G. Kanbach for making us aware of this possibility). 

In the presence of prolonged capture, the DS could end 
up even larger, possibly leading to supermassive stars 
that consume all the baryons in the original halo, i.e., 
possibly up to 10^ Mq stars. The role of angular momen- 
tum during the accretion must be accounted for prop- 
erly and will affect the final result; in fact, even for the 
IOOOM0 case studied in this paper, we intend to examine 
the possible effects of angular momentum. This possibil- 
ity will be further investigated in a future publication (in 
preparation), both from the point of view of the stellar 
structure, as well as the consequences for the universe. 
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